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sis b ten?orhit^r k u ^ S of tbe m,ddIe ma gnetotail were developed using a con- 

orrlnnc C ? techn,t l ue ; II was found that currents carried near the equator by 

groups of ions with anisotropic distribution functions are not well approximated by the 

fated" 8 center ex P ress ' on s. The guiding center equations fail primarily because the calcu- 

tne P r SUre t6nSOr ,S n0t T agnCtiC fi6ld aHgned - 7,16 P ressu " e tensor becomes field 
a gned as one moves away from the equator, but here there is a small region in which the 

die nrel Ce " r eqUatl ° nS madequate because tbe two perpendicular components of 
prwesses [Lnwa f T UnCq “ a, n The s, S" ificance of nonguiding center motion to substorm 
sheeMnvolves field , mechanism that may disrupt a thin cross-tail current 

est altimd i! fi ' d hangCS u hat cause ,ons t0 be S in following chaotic orbits. The low- 
fmnort!n? TB ha0 reg '° n ’ charactenzed b y an adiabaticity parameter k = 0.8, is especially 
Z e , tf aVerage cr0 f: tai1 partic,e drift is s,ow - and we were unable to generate a 

sheetfiTT ShCe “T 8 SUCh ’T Therefore an y P rocess that tends to create a thin current 
sheet in a region with k approaching 0.8 may cause the cross-tail current to get so low that 

10 SUP? °u the u' 0beS - A diffCrent ,imit may ** ^ant ^ resonant 
orbU^ofath^ current sheet because particles reach a maximum cross-tail drift 

thfSml 1* "fT ber 0f r per un,t length decreases as the tail is stretched, this part of 
the ln£T TV 6 if ? may U " able t0 Carry the cross ' tail current needed to support 

the dt5h T f hin f Hee afe neCded for both resonant and chaotic orbit mechanisms because 
the distribution function must be highly structured. A description of current continuity is 

included to show how field aligned currents can evolve during the transition from a y 

two-dimensional (2-D) to a 3-D configuration. 


L Introduction 

A companion paper [Kaufmann et at., this issue] and Larson 
and Kaufmann [ 1996] (hereafter called LK96) describe the consis- 
tent orbit tracing (COT) technique used here. The procedure starts 
by picking magnetic and electric field models. Orbits of many 
groups of particles are traced in these fields to calculate the electric 
current each group carries. The groups then are combined so that 
ions and electrons in the full plasma carry the electric current 
needed to generate Ihe preselected magnetic field. 

The region studied is -20 R E < x < - 14 f? £ , 0 < Izl < 2 R F in geo- 
centric solar magnetospheric (GSM) coordinates. The field model 
is symmetric about the equatorial plane, z = 0. The standard mag- 
netic field from LK96 was used for most of the calculations in this 
paper. The model suggested by Him et at. [1975] and Zwingmann 
11983] provided the dominant contribution to our standard model 
in the region of interest. Parameters were adjusted to give a charac- 
teristic current sheet thickness of about 0.7 R f and so that both the 
total cross-tail sheet current density and the adiabaticity parameter 
would be similar to those obtained in the Kp = 4 version of the 
Tsyganenko [1989] (T89) model. The adiabaticity parameter is 
defined by k = /? min /p max [Buchner and Zelenyi, 1 989], where 
^min | s minimum magnetic field line radius of curvature and 
Pmax is ,he maximum particle gyroradius. Both R and o 

min rmax 
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noo f f Und at 2 = °- A di P° le and lhe Tsyganenko and Usmanov 
[ 982] ring current model produced the dominant magnetic fields 
at lower altitudes. 

The output from the COT analysis is a three-dimensional (3-D) 
velocity distribution function /( r, v) for each spatial box. The 
boxes are Az = 0.03 to 0.1 R £ wide in the z direction and Av = 0.5 
or I R e wide in the a: direction. Only one wide box was used in the 
y direction, and the model is referred to as two dimensional 
because there is very little y dependence in the region of interest. 
The calculated /( r, v) then were used to evaluate fluid parameters 

such as the mass density p, the bulk velocity V, and the pressure 
tensor P. 

Ion orbits in the neutral or inner plasma sheet region deviate 
markedly from the spiral guiding center motion that is typical of 
the outer plasma sheet and the radiation belts. Nonguiding center 
orbits have been studied by many groups. Previous studies were 
summarized by Kaufmann et at. [this issue] and by LK96, so we 
will refer to these papers rather than repeating the summary. 

The goal of section 2 is to find which bulk parameters are 
strongly dependent on nonguiding center properties of particle 
orbits in the middle magnetotail. Ft was shown why the full pres- 
sure tensor is needed near z = 0 rather than only the pressures par- 
allel and perpendicular to the magnetic field. The calculated 
pressure tensor is not oriented along B in the weak field region near 
the equator. A region also exists in the model in which P is nearly 
field aligned but the perpendicular pressure is not the same in all 
azimuthal directions around a field line. 

In section 3 the results from the nonguiding center study are 
applied to substorm problems. It is concluded that cross-tail current 
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may be disrupted in a region which is initially dominated by 
Speiser particles [ Speiser , 1965]. Speiser particles follow spiral 
guiding center orbits while they are well away from z - . ese 
particles carry cross-tail current very efficiently at Izl < z 0 where 
they meander one or more times across the equatorial plane. The 
location z = z is the point at which the z component of the parti- 
cle’s gyroradius is equal to the particle’s distance from the equato- 
rial plane. This thickness parameter is defined by z 0 
mv/ lqB„ (z o) ] , where s] y = B ] + B 2 y ,m is the particle mass, v 

is the velocity, and q is the charge. 

One specific disruption process involves a thin current sheet 
region in which the thickness and/or B. changes until most current 
carriers have an adiabaticity parameter near K = 0.8. Most particles 
with ic near 0.8 suffer rapid chaotic changes in orbital properties 
each time they interact with the neutral sheet, so the ions o no 
follow Speiser orbits for even a single interaction. We were unable 
to generate a consistent thin current sheet using ions wit K near 
0.8 The next chaotic region at K - 0.4 was found to support a thin 
current sheet. 

Particles can remain on Speiser orbits and drift rapidly across 
the tail in a thin current sheet when their average k parameter is 
near a resonant value such as 0.5. Here current disruption can be 
produced by a decrease in the number of particles in a resonant 
region because particles already are drifting at their maximum 
cross-tail speed. The fewer particles per unit length of the tail are 
unable to support the lobe field. Finally, reasons for changes in 
plasma sheet thickness and the relationship of 2-D steady stat 
COT results to a 3-D time dependent tail are discussed. 

2. Nonguiding Center Effects 

The relatively thick standard model from LK96 will be used in 
this section. A uniform 0.3 mV/m cross-tail electric field was 
included to produce earthward plasma drift. The COT method 
traces orbits of a group of particles with preselected initial energies 
starting from a point that was randomly selected within one of the 
spatial boxes. Particles are energized by the electric field because 

they also drift in the cross-tail direction. 

The inner current sheet is of particular interest because the 
model magnetic field is too weak at Izl < z 0 to deflect ions through 
90° back toward the equatorial plane. This section wi sow 
the pressure tensor is oriented in directions that are unrelated to the 
magnetic field and that the strength of the cross-tail current devi- 
ates substantially from guiding center predictions in the inner cur- 
rent sheet. In the outer current sheet beyond about Izl - 2z 0 , ions 
spiral around the magnetic field line direction, and the guiding cen- 
ter approximations are valid. 

2.1. Nonguiding Center Particles 

LK96 showed that nonguiding center magnetotail ion orbits can 
be divided into three classes according to the cross-tail current dis- 
tribution that each carries. This separation depends on the points at 
which particles either mirror or are reflected back toward the equa- 
tor. Figure la is the y-z projection of a small segment of a particu- 
larly simple nonguiding center trapped particle orbit An ion that 
follows this figure eight shaped trajectory as it meanders ac an 
forth across the equatorial plane carries positive j y at z > z Q an 
negative j at Izl < z„. There is only a small net drift in the positive 
y direction for the full orbit shown yielding a small average drift 
velocity v . The subscript g indicates that this orbit-averaged 
drift corresponds to the bounce-averaged drift velocity of guiding 
center particles. The large localized currents associated with 


motion in the positive direction at Izl > z Q and negative direction at 
Izl < z correspond to magnetization currents j my = nqv my , where n 
is the number density and q is the charge. The total cross-tail veloc- 
ity is the sum of the magnetization and orbit-averaged velocities, 


1 A + second class of particles mirrors in the outer current sheet, 
beyond Izl = 2z_. These trajectories were called cucumber orbits 
by Buchner and Zelenyi [1989]. Figure 2 of LK96 shows particles 
following cucumber orbits. Finally, the ion in Figure lb which mir- 
rors far from z = 0 is on a Speiser orbit. Such particles carry most 
of their cross-tail current in the +y direction at Izl < z 0 and produce 
the thinnest current sheets. 

The solid curves in Figures 2a and 2b show the current distribu- 
tion j v (jr, z) carried in the region -17 R £ < x < -15 0 < Izl < 

R by a group of 1 000 protons that were randomly selected from an 
isotropic 5-keV Maxwellian parent distribution starting at U, y, z) = 
(_l 55 ^ o, 0) R . Figures 2a and 2b show the current distribution as 
a function of z within a Ax = l-/? E -wide a- box. Although few of 
these 1000 trajectories were as simple as the one shown in Figure 
la, the group was dominated by particles that were trapped at Izl < 
2z =05 R p while in these x boxes. The pattern of negative/^ at Izl 
< z , positive j at Izl > z„, and small j in the outer current sheet is 
typical of trapped particle groups [ Ka'ufmann and Lu , 1993]. Runs 
also were carried out using smaller z boxes to examine the interest- 
ing nonguiding center region in more detail. Figures 2c and 2d are 
these expanded plots for a similar randomly selected group of ions^ 
The solid curves in Figures 3a and 3b show j y (*, z) in the 
/? < * < -17 R , 0 < Izl < 0.8 R e region for a group of 1000 pro- 
tons selected from a 5-keV Maxwellian distribution that started at 
(-9.05, 0, 2.05) R f in the outer plasma sheet. Most of these parti- 
cles follow Speiser orbits. Their typical current pattern is concen- 
trated at Izl < z = 0.25 R e but often exhibits a minimum very close 
to z _ o Finally, the solid curves in Figures 3c and 3d show 
j (x, z) carried in the -19 R E < x < -17 R E , 0 < Izl < 2 R f region 
by ions and electrons in the combined plasma that was generate 
by the COT analysis. This current nearly equals the j y (x,z) 
required to generate the standard model magnetic field. 

2.2. Guiding Center Particles 

The guiding center drift equation for particles of species s is 


i - — x 
J ,c - fi2 * 


VP 


A'l 




+ nA ^B + „ m! B x[( v 5 .V)V s ] 
B B 


( 1 ) 


where parallel and perpendicular components of the current j t 
and pressure tensor P s are taken with respect to the direction 
of the magnetic field B. Equation (1) can be obtained by com- 
bining guiding center gradient, curvature, E x B . polarization, 
and magnetization drift terms [Parker, 1957]. This expression 
was used to calculate electron current in the COT analysis. 
The term in (1) involving the bulk velocity V, is the polariza- 
tion drift and is small in this magnetotail model. 

A generalized version of the steady state drift equation 

i _ H x V dP '- P° 
lsl B 2 (jT, 3 *P 

+ „Ex B B [ (Vj .V)V s ] (2) 

s * B 2 
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Figure 1. (a) A simple nonguiding center trapped particle orbit, (b) a Speiser orbit, (c) a guiding center orbit show- 
ing cross-tail magnetization currents, and (d) a guiding center orbit showing cross-tail drift motion. 


is obtained by taking the cross product of the momentum 
equation [Kaufmann et a/., this issue, equation (6)] with B and 
dividing by . The subscripts a and p refer to Cartesian 
components in GSM coordinates, and the sum is the diver- 
gence of the pressure tensor V • P v . Equation (2) follows 
directly from the Vlasov equation and should be valid for any 
group of particles in such a steady state plasma. Equation (2) 
reduces to (1) when can be written in the form 

P, = /VJ + (/>., n -P lX )bb (3) 

where b is a unit vector along B and 1 is the unit tensor. These 
approximations are valid for electrons but not for ions in the 
inner current sheet. 

Figure Ic shows the x-z projection of a trapped particle that fol- 
lows a guiding center orbit. Cross-tail magnetization current in the 
— y direction near z = 0 is primarily given by the second P s ± term 
in the square brackets in (1) and is physically produced by the 
crowding together of orbit segments on the concave side of the 
field line in Figure 1c. This particle is always moving in the +y 
direction when it is farthest from z = 0, resulting in a positive / . 
Figure Id is the y-z projection of a guiding center particle that mir- 
rors far from the equatorial plane. The cross-tail motion of the par- 
ticle’s guiding center near z = 0 is produced primarily by magnetic 
field curvature drift or the P s ^ term in (1). 

Figures la - Id illustrate some of the physical distinctions and 
some similarities between currents carried by guiding center and 
nonguiding center particles. Although the orbits in Figures la and 
lc look quite different, both involve strong magnetization currents 
flowing in the +y direction at the largest Izl reached by the particles 
and in the -y direction near z - 0. In both cases there is more back 
and forth motion than net drift, indicating that v > . 


One important distinction is that the thickness parameter z Q has 
special significance only for the nonguiding center ions. Many ions 
in the randomly selected groups of 1000 that produced the 
j y (jc, z) distributions in Figures 2a - 2d are on trapped orbits that 
meander back and forth across z = 0 while remaining at Izl < 2z 0 . 
The solid curves in Figures 2e and 2f show currents carried by a 
guiding center ion group. These particles also were started at 
(-15.5, 0, 0) R e with an isotropic distribution. An energy of 10 eV 
was used so the ions would have such a large K that they follow 
guiding center orbits. Although the structure of negative j at z = 0 
and positive j at large Izl is qualitatively similar to that of trapped 
nonguiding center ions, the scale is much different. The z Q parame- 
ter has no significance for guiding center particles because they fol- 
low spiral rather than meandering orbits. The mirror points of 
guiding center particles depend only upon their equatorial pitch 
angles and the magnetic field structure, not upon their energy or the 
energy dependent z Q parameter. 

Figures 2e and 2f show that both the dotted curve calculated 
using (1) and the dashed curve calculated using (2) produce good 
estimates of j ( x , z) for guiding center particles. The structure of 
j y (jc, z) for these nearly guiding center particles was found to 
have a weak energy dependence. Figures 2e and 2f which used 
1 0-eV protons have peaks at z = 0. 1 3 R £ . A similar plot was gener- 
ated using 1-eV protons, and the peak moved to 0. 10 R^. If the def- 
inition of z in section 1 is used, then z = 0.01 1 and 0.0035 R r for 
10-eV and 1-eV protons, respectively. Figures 2e and 2f therefore 
show that z Q is not relevant for guiding center particles. 

Equation (1) shows that j y (x, z ) depends on temperature and 
density or pressure variations for guiding center particles. The den- 
sity and temperature of particles vary as they drift in the electric 
field. However, in a 2-D model the distribution of particle energy 
as a function of jc is independent of the magnitude of £ y , provided 
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Figure 2. A comparison of j determined by integrating /( r, v) with the approximations given by equations (1) 
and (2) for two x boxes. Ions were injected with an isotropic distribution function at ( x , y, z) = (-15.5, 0, 0) R E pro- 
ducing a group with many trapped ions, (a) and (b) The value j carried by 5-keV nonguiding center ions in the 
range 0 < Izl < 2 tf £ . (c) and (d) More detail in the 0 < Izl < 0.6 /? £ region, (e) and (0 are similar plots for 10-eV guid- 
ing center ions. 


E y * 0. For example, a particle that starts with 5 keV of energy at 
the equator at x = -20 fl £ will be more energetic when it has drifted 
earthward and crosses the equator at x = -14 R E in a uniform 
cross-tail field E y = 0. 1 mV/m. The particle gains energy because 
the magnetic field causes duskward drift. Doubling £ y to 0.2 mV/m 
doubles the earthward drift speed but does not change the cross-tail 
drift speed at a given x in the 2-D model. As a result, the particle in 
the larger E y will reach x - -14 R £ twice as quickly, drift half as far 
in the y direction, and therefore gain the same energy in the 0.2 
mV/m electric field as a particle that drifts to the same jc in an 0.1 
mV/m field. Similarly, the ratio of ion to electron energies remains 
constant if they E x B drift earthward and their cross-tail drift 
speed depends only on energy. 

Figures lb and Id show that both guiding center and nonguiding 
center particles with mirror points near the Earth carry their stron- 
gest cross-tail current near z = 0. However, only the nonguiding 


center or Speiser particle current is confined to the region Izl < 2z 0 , 
and only the nonguiding center particle j ( x, z) shows the reso- 
nant structure that will be described later. The solid curves in Fig- 
ures 3e and 3f show the current carried by IG-eV guiding center 
protons that were injected at the same point as the 5-keV nonguid- 
ing center particles in Figures 3a and 3b. The qualitative similarity 
and difference in scales again is evident. The guiding center parti- 
cles all have mirror points earthward of the (-9.05, 0, 2.05) R 
injection point so the distribution function is field aligned at z = 
and F|| > P L . As a result, the second term in (1 ) is dominant in the 
current sheet. 

2.3. Breakdown of the Guiding Center Approximation 

The solid curves in Figures 2, 3a, 3b, 3e, and 3f are nqV for 
single groups of 1000 ions each, where the bulk velocity V is cal- 


Ot*i 
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Figure 3. Similar to Figure 2 except (a) and (b) the 5-keV ions were injected at (jr, y, z) = (-9.056, 0, 2.053) R E pro- 
ducing a group of Speiser particles, (c) and (d) Results from the combined consistent orbit tracing (COT) distribution 
function, (e) and (0 Similar to Figures 3a and 3b except 10-eV guiding center ions were used. 


culated by integrating over the COT distribution functions 

f s {r,v) 

^.v.a = ^ /»’„/, (r, v )^ v (4) 

The solid curves therefore show the actual current carried by 
all particles in the model plasma. The individual-pressure ten- 
sor [Rossi and Olbert , 1 970] 

P l% = m s\ [v 0 - [vp- V, p] 4 (r, Y)d\ (5) 

was also calculated by integrating over / Y (r, v) . The super- 
script in (5) indicates that is calculated in a frame mov- 
ing at the bulk velocity of the species s rather than in the frame 
moving at the bulk speed of the entire plasma. This distinction 


is needed when looking at individual ion groups. Each ion 
group is treated as a separate species so that (4) and (5) give 
the bulk velocity and partial pressure tensor for that group. 
Since (2) requires taking numerical derivatives, the current 
calculated using this expression is less accurate than that 
shown by the solid curves. All of Figures 2 and 3 show good 
agreement between the kinetic results based on (4) and the 
dashed curves which are based on (2). The point at smallest izl 
has been omitted from all plots derived from (1) and (2) 
because the numerical derivatives have larger uncertainties at 
these end points. 

The dotted curves in Figures 2 and 3 were calculated using (I) 
and show how much this guiding center expression deviates from 
the actual current carried by groups of anisotropic ions. These and 
similar figures for other particle groups showed that the three meth- 
ods used to calculate j y (.v, z) agree beyond approximately Izl = 
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2z 0 = 0.5 R e . Causes of the differences at Izl < 2z 0 are described 
below. Figures 3c and 3d show that (1) and (2) give reasonably 
good estimates for j y in the full consistent thick current sheet. This 
agreement is a result of the near isotropy of the combined ion dis- 
tribution function near z = 0 in the thick model current sheet. 

2.4. Pressure Tensor Elements 

In any 1-D magnetotail model with fi y = 0 and 3/ 3x - 3/3y = 

0, the only nonzero contributions to V P are 3P../ 3z and 
3 P t _/3z . The P.. derivative produces particle forces in the z 
direction which balance the magnetic force associated with the 
increase of B as one moves from z - 0 to the lobes. The P xz term 
provides the particle force that balances the magnetic force in the x 
direction. 

Each element of P involves an integral of the distribution func- 
tion times a product of two thermal velocities , where w a = 

v - V a . The bulk velocity V a is small compared to the total par- 
ticle speed v Q and to for most particles in the quiet time 
model. Diagonal terms in the pressure tensor are largest because 

is positive at all velocities. Off-diagonal elements such as P x: 
would be very small if and w, were uncorrelated since velocity 
space regions with positive and negative w x w. products would 
nearly cancel. However, these two components can be highly cor- 
related in a magnetized plasma. 

As one simple example, consider a guiding center particle that 
mirrors near the Earth. When this particle is approaching the equa- 
torial plane from its northern hemisphere mirror point, it will have 
an average w < 0 and < 0. When the particle moves away from 
the equatorial plane toward its northern hemisphere mirror point, it 
will have an average w > 0 and w m > 0. The w x w : product there- 
fore is usually positive near the equator and can produce a substan- 
tial off-diagonal pressure tensor element P x: . The product is 
largest when w and are approximately equal, which corre- 
sponds to the point at which I B\ = \B J for a field aligned guiding 
center particle. The physical significance of a large P xz when 
3/3z » 3/3x or boxes are selected with Ax » Az is that particles 
which move from one z box to the next carry substantial x momen- 
tum. If 3P v _/3z is large, then particles which enter one side of a z 
box and leave the other side will deposit net x momentum in the 
box. For the field-aligned particle example used above, particles in 
the northern hemisphere that are moving toward z = 0 tend to have 
a larger tailwardx velocity when they enter a z box than when they 
leave it. Particles moving away from z = 0 tend to have a larger 
earthward x velocity when they leave a z box than when they enter 
it. In both cases the particle has gained momentum in the positive x 
direction while inside the z box. This requires an earthward force to 
be exerted on the particles. An equivalent physical picture is that a 
tailward centrifugal force is exerted by the particle as its guiding 
center moves along the curved magnetic field line. 

Other pressure tensor elements also contribute to force balance 
in the 2-D model used here. The off-diagonal pressure tensor ele- 
ments can be 10% or more of the diagonal elements for individual 
particle groups. Equation (2) shows that the x and z components of 
V • P contribute to j . The 3P„/dz term generally dominates 
over the 3P r ./3x term in the z component of the divergence, just 
as for the 1-D case. However, dP XK /dx and 3P r ,/ 3z are compa- 
rable in the x component of the divergence in the interesting region 
near z = 0. For studies of j x and j z , it is the 3P yr / 3z term that is 
usually dominant in the 2-D model. 

2.5. Normal Coordinates 

Guiding center particles are most naturally described in a mag- 
netic field oriented coordinate system. For a mirroring particle the 


positive and negative products of the guiding center thermal veloc- 
ities parallel and perpendicular to B are equally likely, so the 
w \\ w i P r °d ucts and associated off-diagonal pressure tensor ele- 
ments average to zero. The thermal velocities w'jj and alcmg 
any two orthogonal directions perpendicular to B are highly corre- 
lated but are 90° out of phase for guiding center particles so the 
average w L] w 12 product and associated off-diagonal term in the 
pressure tensor also is zero. As a result, only the diagonal P n and 
P u = P 12 = P ± terms are nonzero as is assumed in (3). 

Since P is a real symmetric tensor even for nonguiding center 
particles, it can always be diagonalized to give real eigenvalues 
and eigenvectors or normal coordinates. The three nonzero ele- 
ments of P in the nonfield-aligned normal coordinate system will 
be called P, , P n , and P m . There are two ways in which the 
assumptions that give (3) can be violated. First, it is possible that 
none of the normal coordinates of P is along B. Second, even if 
one normal coordinate is nearly along B so that P ] = P\\ , it is not 
necessary for the other two elements P n = P ±] and Pm = P±2 
to be equal. 

Diagonal ization and the calculation of normal coordinates relies 
upon accurate values of the off-diagonal elements. Distribution 
functions from the COT analysis produced the most reliable nor- 
mal coordinates in the inner current sheet. LK96 showed that the 
ion temperature at the outer edge of the model current sheet is 
about half the temperature near z = 0. At large Izl it is often found 
that the calculated velocity distribution function is zero in many of 
the higher velocity boxes because few of the ions which were 
traced had large enough velocities at this location to make contri- 
butions to the high-velocity boxes. Moment equations such as (4) 
and (5) give reliable values of the density, bulk velocity, and the 
large diagonal elements of the pressure tensor throughout the cur- 
rent sheet. However, the small off-diagonal elements of P are less 
reliable at large Izl. 

Figure 4a shows the angle 0 between B and the closest normal 
coordinate axis for trapped ion groups with the same starting 
parameters as those in Figure 2. Figure 4d is a similar plot for 
Speiser groups with the same starting parameters as those in Figure 
3. Five runs using different randomly selected groups of 2000 ions 
each were used to produce the five curves that are superimposed in 
Figures 4a - 4f. This presentation was used to show the reliability 
with which normal axes can be calculated using the available 
/( r, v) . The runs used 2000 ions each to produce more accurate 
off-diagonal elements and Az = 0.05-P^-wide z boxes in the 0 < Izl 
< 1 R e range to emphasize detail in the inner current sheet. 

An inspection of the GSM components of the normal coordinate 
axes derived by diagonalizing P shows that near the equator the 
normal coordinates were consistently near the GSM x, y, and z 
axes. This feature is produced by the symmetry of the magnetic 
field model and of the injection scheme that was used. Since the 
magnetic field is symmetric about z = 0, distribution function infor- 
mation was folded about the equatorial plane. Particles contributed 
to the distribution function in the same spatial box whenever they 
were located at a distance Izl either above or below the equator. 
Physically, this corresponds to starting pairs of particles, one at a 
point (x, y, z) and another at (x, y, -z). 

The magnetic field is in the z direction at z = 0 and points closer 
to the GSM z axis than to the other two GSM axes at Izl < 0.1 
the two smallest z boxes in Figure 4. Since the normal coordinates 
of P at Izl < z Q = 0.25 R £ are nearly along the GSM x, y, and z axes, 
the angle 0 is roughly equal to the angle between B and the z axis 
for the first two points on each curve in Figures 4a and 4d. The 
magnetic field points closer to the x GSM axis in all other z boxes. 
The angle 0 is roughly equal to the angle between B and the x axis 
for the next three smallest z boxes in Figures 4a and 4d. The fifth z 
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the two diagonal pressure tenfor witi^lSSISir”" 5 ‘T™ "T 3 ’ VeClor (b) and (e) The ralio of 

that involves all diagonal pressure tensor elements. This raliJT? /?“ ifrtlrT f ar,hes ‘ 1 from B < c > and (0 A ratio 
The ions in Figures 4a - 4c were injected at (-\ <5 S n m y? J 1 • 1 ^ first normal coordinate points along B. 

4d- 4f were injected at (-9.056, 0, 2.053) R £ producing a gro^ofS^is^Hons.^ ‘ raPPCd ^ ^ ' h ° Se m FigUreS 


™tes nf P V 0 ' y ° nd ,his P° int one of the no mal coordi- 
nates of P points closer to B than to any of the GSM axes so the 

pressure tensor has become field aligned. Figures 4a and 4d show 
that 9 becomes very small, indicating that one of the normal coor- 
dtnates is very close to B. somewhere in the range z Q < |-| < 2 z 0 = 

To be more quantitative, recall that the 2-D model magnetic 
field .s m the x-z plane. The x-z plane projection of one normal 
oordmate was found to point within about 1° of the magnetic field 
direction everywhere beyond Izl = 1 .5 z Q for the Speiser group. This 
normal coordinate fluctuates by 2° to 3° out of the x-z plane P How- 
cverthesesmalldeviations in the y direction appear to be edge 

feren T * Par ' iCleS e " tering and leavi "g * box at dff- 

ferent z locations. The y deviations become smaller as particle 

efferr" 8 p”"' 8 “? Spread over larger regions t0 smooth out edge 
effects. Figure 4a shows that the trapped particle group also 

becomes field aligned near Izl = 1.5 z o . The increasing scaUer in 
gure 4a near Izl = I R £ shows that the pressure tensor elements 
are not accurate enough in this region to calculate reliable normal 
coordinates for the trapped particle group. The trapped parole 
results are inaccurate near and beyond I R because few particles 
m the group reach this far from z = 0. 

haveteen or/"^ ^ ^ P " /P ™ where the e,e ™nts 

have been ordered so that P„ > , n the region Izl > 1.5 z 

is eld aligned, P u and P f[| are the two perpendicular 

pressure elements. Comparing Figure 4a with Figure 4b and Figure 

4d with Figure 4e shows that there is a region somewhere between 

V?.?° m Wh,Ch P iS field aligned or 0 is near| y zero but in 

which there is substantial perpendicular asymmetry (/>., * P ) 

S t “ “ ‘ ^ “ »'« ■'*»! a"- "ymmSic 

" H ” F ; g, "° S « *" d *= tot the average calculated 

nni'n. I! 'h 1,0 devia,es from one by no more than the scatter of ' 

d P Zt ' C r g thiS devia,i0n can be ex P ,ained by errors in i 
determining the normal coordinates. 

Figures 4c and 4f show 2 />,/(/>„ + />,„) which is P„/P, 


en P is field aligned and symmetric about B. The qualitative 
behavior of the pressure ratios agrees with expectations for the 
orbit groups selected. Particles should have P,, <p near their 

P ° in,S - M ° SI Par,idcS in the general| y mapped 
g p -igure 4c) are reflected in the region shown. The Speiser 
particles mirror beyond the region plotted in Figure 4 Their 

P " /P \ ratl ° (Flgure 4 0 is expected to increase through the ordi- 
nary mirror effect as particles move toward the equator as long as 
the magnetic moment p is conserved. The Speiser particles become 
more nearly isotropic a. the smallest Izl where p is no. conserved. 
Although qualitative features can be understood, a quantitative 
analysis is difficult using the information in Figure 4. It must be 
recalled that each panel shows variations as r increases at a fixed v 
not the variation seen when one moves along a field line. Particles 
cate al - I R ^ in the thick model cross the equator 3 to 4 R 
tail ward of the ,v box being studied. " E 

To summarize this section, it was shown that the direction of B 
is almost irrelevant to the orientation of P at Izl < z for groups that 
are dominated either by trapped or by Speiser particles Section 3 
descntes a thin current sheet that was made using only Speiser par- 
ticles. It therefore is not sufficient to calculate only P., and P in 
a thin current sheet when analyzing satellite measurements near z = 

0 All I elements of P are needed. For the cases studied, the calcu- 
lated P consistently became field aligned somewhere in the region 
‘o < - 1 < 2 V A small region was found extending beyond the 

P0 ' n ; a ' wh ' Ch P became ,ield aligned in which (I) and (3) are 
invalid because P„ * />,„ (Figures 4b and 4e). I, was found that 

' , t (3) c . at1 bc u " cd be y° nd W = 2z 0 . These conclusions sug- 

gest hat evaluating P may provide a way to determine whether* 
satellite is near z - 0 in a thin current sheet. If the observed P is not 
field aligned then the satellite is likely to be located at Izl < 1.5 z 
he test is not useful when B is nearly along a GSM axis since then 

' “ . n °‘ P° ss,b *e to distinguish between GSM-aligned and 

field-aligned orientations. 

Figure 5 shows the same parameters that were plotted in Figure 
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Figure 5. Si™,., » Figura 4 — 

(c) A ratio iha, involve. .11 diagonal prae.ura ienso, elements. This ra„o 
r, / r. ■£■ rr\nrdinate noints fllong B. 


axes maL die launvov • vy 

P\\/P 1 if the first normal coordinate points along B 

4 but now for the final /( r, v) that was produced by summmg ion 
groups to produce the thick COT model. Figure 5a shows that 10 
becomes small or P becomes field aligned near Izl - 1.3 z ■ The fu 
0 < Izl < 2 R F range is plotted in Figure 5 to show that the normal 
coordinate determination fails so the calculated Qbecomes erratic 
near |z | - 1.5 R The z boxes are twice as large in Figure 5 as in 
Figure 4. As a result, the 6 variations very close to z = 0 are not 
well illustrated in Figure 5. However, the normal coordmates again 
ate aligned with the GSM axes near z = 0. Figures 5b and 5c show 
that P is nearly isotropic and nearly symmetric in the region Izl < z Q 
since P, . P„ . and P,„ ara all hearty equal. A number o auihors 
have shown that P must approach isotropy near Izl = 0 if /(r, v) 
symmetric in an equilibrium current sheet [Cole a c l " ^ 
1972; Rich el al„ 1972; Cowley, 1978; Nonet el at., 1985, Hill and 

Voigt, 1992 ], 


3. Substorm Effects 
3.1. Kappa Dependence 

The observation that thin current sheets often form before sub- 
storm onset suggests that nonguiding center effects may play a role 
TZ substorm process. Orbits were classified in section 2 accord- 
ing^ particle minor or reflection points because this separation 
info trapped, cucumber, and Speiser orbits determines the basic 
shape of the cross-tail current distribution j (x, z) . Processes di - 
£Ert iu rt»> section also depend nn ihe dynamical ebar.c 
Of particle orbits. U is possible for a group of Speiser particles to be 
dominated by either resonant or chaotic orbits. The omin 
dynamical characteristic of a group of particles depends on the K 
parameter (section 1), which varies with both the ion energy and 
fhe current sheet thickness. The first and second orbital resonances 



ages for particles that cross z = 0 at this * location. 
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appear near k = 0.5 and k = 0.3, respectively. Most particles in a 
group with K = 0.5 cross the equatorial plane 2 times (Figure lb), 
and most particles with k = 0.3 cross the equator 3 times per cur- 
rent sheet interaction. Ions in a group with a resonant k parameter 
leave the current sheet region with almost the same magnetic 
moments they had when they entered. Almost all resonant ions that 
enter the current sheet on Speiser orbits will also leave the current 
sheet on Speiser orbits. When k is near 0.4 or 0.8, particles exhibit 
their most chaotic behavior. The mirror points of chaotic particles 
that leave the current sheet depend very sensitively upon all orbital 
parameters of the particles when they entered. An incident group of 
chaotic Speiser particles becomes thoroughly mixed between 
Speiser, trapped, and cucumber orbits. 

Section 2 was based on studies of a thick quiet time current 
sheet. The characteristic magnetic field z scale length L w was much 
larger than the thickness parameter z o (section I) for 5-keV pro- 
tons. Individual groups of particles had highly anisotropic distribu- 
tion functions near z = 0, but the final consistent current sheet was 
nearly isotropic because the consistent thick model required a mix- 
ture of all orbit types. This section emphasizes results from a 2-D 
equilibrium model of a thin current sheet (A. Bhattacharjee, per- 
sonal communication, 1996). A current sheet with thickness L^z 
is important because it can be created using only Speiser particles. 
The resulting consistent thin current sheet has a highly anisotropic 
distribution function. 

Figure 6 shows the most significant result obtained from study- 
ing thin current sheets. The plot shows contours of the parameter 
VV w here V is the cross-tail ion bulk velocity and v is the 
total ion velocity. Groups of monoenergetic 5-keV Speiser protons 
were started at z = 2.5 R £ along a number of field lines to prepare 
this plot. Monoenergetic ions were used because ic is energy depen- 
dent, and we wanted all ions in a group to have nearly the same k. 
Ions drifted through the region of interest, but V was retained only 
in the one or two 0.5-tf £ -wide x boxes nearest the point at which 
the ions first crossed z = 0. These runs used Az = 0.05-/? £ -wide z 
boxes, and 15 ion groups were traced to create Figure 6. 

Note that V y /v 0 in the -16.0 R £ < x < -15.5 R f x box exceeds 
0.3 at z = 0, drops to 0.2 at z = 0. 1 /? £ , and drops to zero at z = 0. 1 5 
R f . This concentration of current into a well-defined sheet at Izl < 
z o = °* 1 R e is l yP ical ° r Speiser orbits (Figures 1 b, 3a, and 3b). The 
K parameter defined in section 1 was determined by calculating the 
field line radius of curvature R /R 1 = -(B V) B /B 2 and the 
maximum ion gyroradius p max = mv 0 /qB each time an ion 
crossed z = 0. The K values shown at the top of Figure 6 are aver- 
ages for all particles that crossed the equatorial plane at a given x 
and therefore refer to ions near z = 0. A particle that is located at z 
= 1 /? £ in the -16.0 /? £ < x < -15.5 /? £ x box will cross the equato- 
rial plane 3 R £ tailward of this region and will have a different K 
(Figure 6). The reason V y /v 0 is so large near z = 0 in this particu- 
lar x box is that K = 0.5 is a resonant value. Almost all these ions 
remain on Speiser orbits throughout their current sheet interaction. 
The maximum value of V y /v 0 for any one Speiser particle would 
be produced if the particle followed a semicircular path in the x-y 
plane with no bounce motion in the z direction. Such a particle 
would have V y /v 0 =0.6. 

The -13.5 R E < x < -12.5 /? £ x boxes in Figure 6 show the 
opposite extreme. Particles crossing z = 0 in this region have an 
average k = 0.8, which corresponds to the most chaotic behavior. 
Although they started at z = 2.5 /? £ , the trajectories do not resemble 
typical Speiser orbits throughout even one current sheet interac- 
tion. The V y /v 0 parameter is less than 0.03 at z = 0, and Figure 6 
shows no evidence of a thin current sheet at this location. 


The -18.0 R e <x< -17.5 R £ x box contains ions with an aver- 
age k = 0.4, which is the second chaotic region. Here V y /v 0 at z = 
0 drops to 0.07, and this parameter reaches zero slightly beyond z = 
0.1 R e . Even though V y /v 0 is small when k = 0.4, Figure 6 shows 
current concentrated in a sheet at Izl < z Q . The second resonance 
takes place near k = 0.3, which occurs in the -20 R £ < x < - 1 9 R £ 
x boxes. Here the V y /v Q parameter again exceeds 0.3, as it did at 
the first or k = 0.5 resonance, and a thin very well defined current 
sheet is evident. Finally, the -1 1 R f < x < -10 R £ boxes show that 
thin current sheets can be produced by the large k nearly guiding 
center particles. 

The monoenergetic groups that were used to prepare Figure 6 
were combined in a COT analysis. Groups with Maxwellian energy 
distributions also were run using the same z = 2.5 /? £ starting 
points. Particles in the Maxwellian groups have a fairly wide distri- 
bution of k values at any one location, producing a mixture of reso- 
nant and chaotic particles. It was easy to find either monoenergetic 
or Maxwellian groups that produced consistent current sheets near 
the resonant regions at k = 0.3 and 0.5. Good models also could be 
produced near the k = 0.4 region of chaotic orbits. When Max- 
wellian groups were used, the k = 0.4 fits were almost as good as 
fits near the resonances. However, we were unable to create any- 
thing resembling a thin current sheet using either monoenergetic or 
Maxwellian particles with an average k = 0.8. Burkhart et al. 
[1992] first noted this problem when injecting particles at large z in 
a 1-D tail model. The significance of this result is discussed below. 

3.2. Disruption of Thin Current Sheets 

Figure 6 shows that any process which uniformly changes k will 
increase the average of nonguiding center ions in one part of a 
thin current sheet and simultaneously decrease V somewhere else. 

A change in either the plasma sheet thickness or in B , changes k 
everywhere. Such a process may induce current disruption if k 
approaches 0.8 in an initially thin region of the current sheet 
because particles with k near 0.8 carry so little cross-tail current. A 
region that is initially thin is needed to cause disruption or current 
diversion by this nonguiding center mechanism because only thin 
sheets are composed of Speiser particles with a highly anisotropic 
/(r, v) . Trapped particles also have a highly anisotropic /( r, v) , 
but current sheets cannot be constructed using only trapped ions 
because the sign of the associated j (x, z) reverses at Izl = z . The 
V v o variati °ns are much less striking than those shown in Fig- 
ure 6, and /( r, v) is much more nearly isotropic when particles 
with different orbit types are mixed as in the thick current sheet 
described previously. 

Since k » 1 in the radiation belts and k generally becomes 
small somewhere in the middle magnetotail [Pulkkinen et al ., 
1992; Kaufmarm et al ., 1993], there nearly always must be a region 
with k - 0.8. Figure 6 shows weak cross-tail current in the +y 
direction throughout a thick region when k = 0.8 even though we 
were trying to generate a thin sheet. This illustrates the tendency of 
particles in a K = 0.8 region to broaden any initially thin current 
sheet. We were able to generate reasonably good COT models of a 
thick current sheet with k = 0.8, though the fits were more jagged 
than the fits with smaller ic shown by LK96 and Kaufmarm et al 
[this issue]. 

A problem arises because a certain sheet current density 
K y (x) , which is j y (x, z) integrated through the thickness of the 
current sheet, is needed according to Ampere s law to produce the 
lobe field. A specific equilibrium lobe field is in turn required so 
that the outward magnetic force per unit area exerted by the lobe on 
the magnetopause balances the total inward force per unit area 
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exerted by the magnetosheath plasma and fields. Kaufmann et al. 
[this issue] concluded that when a current diversion loop forms, 
portions of the plasma sheet outside the loop will be accelerated in 
the tailward direction. The tailward force is generated because 0, 
decreases everywhere outside an intensifying current diversion 
loop. This decrease in can reduce the earthward electromagnetic 
force j y B- so that it becomes smaller than the tailward pressure 
force V P . This will cause a rapid stretching and thinning of por- 
tions of the current sheet [ Ohtam et al., 1992] and an associated 
decrease in k. There may not be enough particles in the thinned 
plasma sheet to carry the required K y if k approaches 0.8, and 
therefore V becomes very small. This process could produce either 
secondary current disruption events or a spreading of the disruption 
region. The initial disruption at substorm onset is discussed below. 

A second important feature of nonguiding center particles that is 
shown in Figures lb and 6 is that the average cross-tail drift speed 
of Speiser particles near z = 0 in resonant regions is a large fraction 
of the total particle speed. This also produces effects that are differ- 
ent from those that would be expected from guiding center consid- 
erations and could be associated with some substorm expansions. 
The guiding center curvature drift term in (l) predicts that j y will 
increase if the current sheet becomes thinner and Pjj > P ± . It 
therefore is possible for the fewer guiding center particles per unit 
x length in a thinning current sheet to continue to carry the K y 
needed to produce the lobe magnetic field. In contrast, Figure 6 
shows that nonguiding center Speiser ions near a resonant region of 
a thin current sheet already are drifting in the y direction as fast as 
is practical. A reduced magnetic field line radius of curvature 
therefore cannot increase V y much, and the change in k is likely to 
reduce V y in such a region. This fact that V y /v 0 cannot be 
increased significantly when it is initially large presents a problem 
if a thin Speiser current sheet is stretched. The tail particles may 
not be able to carry the current required to create the lobe field that 
is needed to balance the solar wind normal force. This process will 
be most important if k starts with an initial value of 0.5 and can 
also be important at higher resonances such as K = 0.3 in the more 
distant tail. It is possible that a sheet can continue thinning until 
anisotropic guiding center electrons carry enough current. Other 
possibilities are that the tail shape may change to reduce the local 
external normal force or the tail could collapse to create a sub- 
storm. 

A final mechanism that can reduce K y without decreasing the 
total number of ions in the tail will be called particle replacement. 
Field-aligned current between the magnetosphere and the iono- 
sphere is a consequence of the need for charge neutrality. Particle 
replacement is a mechanism that can reduce the electron contribu- 
tion to K . The process involves precipitation of energetic electrons 
into the ionosphere and their replacement by low-energy iono- 
spheric electrons. The total number of ions and electrons in the tail 
remains fixed while the K y carried by electrons is reduced because 
low-energy electrons drift slowly in the cross-tail direction. This 
mechanism is unlikely to be important in a thick current sheet 
because ions carry most of the cross-tail current. However, electron 
current could be important on thin sharply curved field lines 
[Pritchett and Coroniti, 1995; Ma et al. , 1995; Birn et ai, 1996]. 
The replacement mechanism may be part of a positive feedback 
effect with precipitation leading to collapse leading to more precip- 
itation. A depletion of the energetic component of magnetotail 
electrons in a restricted region also could be related to the fading of 
arcs before substorm onset [Pellinen and Heikkila, 1978]. 

3,3. Causes of Plasma Sheet Thinning 

Processes that could produce initial current sheet thinning and 
therefore a substorm onset by the nonguiding center mechanism 


are changes in the convection pattern and changes in lobe pressure. 

It has long been known that the current sheet thickness varies from 
several Earth radii to about 0.1 R E while the plasma density 
changes by a smaller factor. Therefore, although some compression 
by enhanced z forces is likely and will be considered below, most 
of the growth phase thinning must involve a removal of plasma 
from the region of interest [Hones et ai , 1971]. The remaining par- 
ticles then must drift faster across the tail to carry the necessary K y . 

Plasma removal implies an imbalance between sources and 
sinks. The primary source of plasma sheet particles for the region 
of interest is inward flow from beyond * = -20 R E . Most of these 
particles are likely to start in the magnetosheath and pass through 
the mantle. Additional contributions come from the flanks and the 
ionosphere [Peroomian and Ashour-Ahdalla , 1995]. Plasma is lost 
from the region of interest through convection toward the Earth 
and eventually out the dayside magnetopause, through losses at the 
flanks, through precipitation into the ionosphere, and through drift 
to open lobe field lines. Since the drift of chaotic particles is differ- 
ent from that of guiding center particles, chaotic orbit effects may 
be able to enhance this drift to open field lines. 

Changes in the interplanetary magnetic field (IMF) are largely 
responsible for changes in the overall growth phase electric field 
and associated plasma sheet convection pattern. It is these 
large-scale convection changes that are likely to be the principal 
cause of thinning the plasma sheet from several Earth radii down to 
< 0.5 R e . For example, the plasma sheet will thin if earthward flow 
continues at the x = -14 R £ boundary of our region of interest and 
decreases at the x = —20 R E boundary. Many substorms are trig- 
gered by the type of changes in the IMF direction which are likely 
to produce changes in the convection pattern [Lyons, 1996]. 

A change in the total lobe pressure near x = -10 R £ also could 
produce a small change in current sheet thickness in the region of 
substorm onset. Changes in the solar wind normal force on the 
magnetopause produce changes in the lobe pressure. Such changes 
are transmitted at the fast magnetosonic speed to the current sheet 
near midnight, where onsets are most common. It is this rapid 
response of the plasma sheet to changes in the solar wind normal 
force that makes this mechanism interesting. The external force on 
the magnetopause is determined by the momentum density in the 
solar wind, by the shape of the magnetopause, and by magneto- 
sheath thermal and magnetic pressures. For example, the magneto- 
pause shape and therefore the normal force are likely to change if 
the location of dayside merging moves from the subsolar region to 
a point near the cusps. 

The observed increase of K y in the inner tail before substorms 
should be associated with some increase in lobe field strength near 
jc = -10 R e . A very crude estimate of the expected change in the 
lobe field is (\i o AK/ 2) (L x /R { ) , where AK/2 Amperes 

per meter of additional sheet current closes in the magnetopause of 
each lobe, L is the dimension in the a direction of the region of 
increased K \ and R f is half the thickness of each lobe. This esti- 
mate assumes that the added current is localized in a region with L x 
< R r Both B x and the lobe flux nB x Rf should increase. An 
increase in flux will cause more magnetopause flaring and there- 
fore more normal solar wind force. Caan et at. [ 1 97 5 1 found such 
changes in lobe pressure during a substorm cycle but did not plot 
absolute field magnitudes. Huang et al. [1992] and Kistler et al. 
[1993] used total plasma plus field pressure measured in the plasma 
sheet to estimate the lobe pressure. These latter studies did not find 
consistent pressure changes during a substomn cycle. 

Although changes in lobe pressure which are produced by sud- 
den changes in the solar wind or IMF cannot produce the large 
changes seen in the thickness of the plasma sheet, they can cause 
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abrupt thickness perturbations. For example, a 20% increase in the 
lobe pressure will compress the high-p plasma in the current sheet 
by about 20%, causing an approximately 10% change in K for all 
particles. The change in current sheet thickness and the change in k 
could be important if the dominant current carrier orbits are near 
either the resonant or chaotic limit discussed above. 

3.4. 3-D and Time Dependent Processes 

Substorms and the diversion of current to the ionosphere are 
3-D, time dependent processes. A 2-D model may be appropriate if 
cross-tail current is disrupted across the tail from flank to flank in a 
limited x region. However, observations show that disruption usu- 
ally is limited in both the x and y directions, producing 
field-aligned currents and a substorm current diversion loop. The 
COT solutions can be used to study some substorm effects even 
when the calculations are 2-D and steady state. For example, forces 
calculated in a 2-D model can provide an estimate of the accelera- 
tion expected near the center of a 3-D current diversion loop. The 
purpose of this section is to examine current continuity as the tail 
evolves from a 2-D steady state configuration to a 3-D steady state 
configuration. The final state includes a confined region of reduced 
V y /v Q . Small steady field-aligned electron currents also are 
present in the final configuration. 

The t < t illustration in Figure 7 shows the initial steady state 
cross-tail current density j y = nqV iy , the ion bulk drift velocity V. y , 
and the ion or electron number density n. The magnitudes of these 
quantities for one fixed x box are plotted along the vertical axis 
while the horizontal axis is the cross-tail direction with positive y 
(dusk) to the left. The three plasma parameters have no y depen- 
dence at t < t Q . Ions drift to the left and electron drift is neglected in 
all but the t > r 4 illustration in Figure 7. No cross-tail electric field 
or earthward drift is included. It is assumed that a region of 
reduced F y /v 0 forms instantaneously over a restricted y region at 
t because something changes K everywhere. The density does not 
change abruptly at t Q , but ions suddenly begin to drift more slowly 
and therefore carry less cross-tail current in this localized region. 
The region of reduced V- is stationary in all the illustrations, and 
various structures drift with respect to this fixed region. 

Figure 7 is not intended to represent the very complex substorm 
expansion process. If j ( x 9 z) suddenly decreased in a large 
enough region of the tail, then a substorm could begin as will be 
discussed in section 4. Here we assume there is only a small 
decrease in V jy and j y in a limited region of the tail. In this model 
the ion density builds up enough by t = t 4 so that j = nqV iy is the 
same as the j that was present at t < t Q . 

Any steady state model requires that electron and ion currents 
must separately be continuous. A transient process involving an 
increasing particle density requires equal ion and electron currents 
flowing into a unit volume to maintain charge neutrality [Atkinson, 
1984]. Since ion cross-tail current drops discontinuously at / in a 
limited y region and we are neglecting cross-tail electron current, a 
field-aligned electron current must be set up during the transient 
phase immediately after r . 

The /, illustration of Figure 7 shows the location of the strong 
field-aligned electron currents that will appear in one to several 
Alfv6n wave bounce periods to close the circuit caused by the dis- 
continuity in cross-tail ion current. This illustration shows the den- 
sity perturbation that is produced because more ions are entering 
than leaving the region with a K associated with chaotic orbits and 
a small V jy . Ions continue to drift rapidly in the +y direction until 
they enter the chaotic K region and resume the rapid drift after they 
leave the region (Figure 6). The density will increase just enough at 
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Figure 7. Sketch of a possible time evolution from a steady state 
two-dimensional (2-D) structure with no y dependence at t < t o to a 
steady state 3-D structure at / > t^. The latter structure contains a y 
region with reduced cross-tail drift speed V jy to show how a con- 
fined chaotic k region could be included in a 3-D COT model. 
Only ion current is assumed to be significant in the current density 
j y , and n is the electron and ion number density. The heavy arrows 
show the strong field-aligned electron flow that is needed during 
the transition between steady states. The light arrows show the 
weak field-aligned electron flow that will be present at the edges of 
the region of reduced V jy for all t> t Q if the slow electron cross-tail 
drift is included. 
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the right side of the region of reduced V jy so that j y = nqV jy is the 
same as at / < / . Ion and electron currents are not separately con- 
tinuous during this transient stage when the density is changing. As 
noted above, there is a net cross-tail drift of ions and a net 
field-aligned drift of electrons into the region of increasing density 
and out of the region of decreasing density. The heavy arrows in 
this and following illustrations show the motion of electrons to and 
from the ionosphere that is required to maintain neutrality. Upgo- 
ing ionospheric electrons and their downgoing current connect to 
the left edge of the region of higher density as this region expands 
to the left. The total number of ions and electrons in the tail is con- 
stant since the density increases in one y region and decreases in 
another. 

The t 2 illustration shows the situation just before the entire 
reduced V y region has become filled with the higher-density 
plasma. The t ^ illustration is still later when the transient region has 
drifted out of the reduced V- region. Eventually, the transient 
region will drift out the dusk flank leaving a new steady state con- 
figuration at t > t A . It is only at this time that there will be any 
change in the total number of ions and electrons in the tail. 

If the slow dawn ward drift of electrons is included, then the 
steady state field aligned electron flow indicated by the light arrows 
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in (he t > f 4 illustration also will be present. The electrons follow 
guiding center orbits so electron drift V does not change abruptly 
at the point the ion k reaches a chaotic value. When cross-tail elec- 
tron drift is included, these small field-aligned currents will be 
present at the edges of the reduced V. region in all but the t Q panel. 
The weak field-aligned current is closed by cross-tail electron drift 
in the higher-density region. It is only the 2-D and 3-D steady state 
configurations t < t Q and t > t 4 that can be modeled using the COT 
method. 

4. Discussion and Summary 

4.1. Nonguiding Center Effects 

The currents carried by groups of nonguiding center particles 
were compared to those carried by similarly injected guiding center 
particles. It was seen that currents carried by the two groups were 
qualitatively similar but that resonant effects and an energy depen- 
dent thickness parameter z Q were important only for nonguiding 
center particles. The model magnetic field in the x-y plane is too 
weak to deflect a nonguiding center particle so that it meanders 
back and forth across the equatorial plane until it reaches Izl = z Q . 
The 2-D field used here had # v = 0. The z Q parameter is smaller 
when a large is added, producing a thinner region of meandering 
particles. 

The most important nonguiding center effects found from the 
COT analysis involved cross-tail currents in plasmas with anisotro- 
pic distribution functions. The frequently used guiding center cur- 
rent equation (1), which is valid when the pressure tensor is 
described by (3), was found to be inadequate in the inner current 
sheet. In this region (2), which was derived from the momentum 
equation and requires knowledge of the full pressure tensor, must 
be used rather than ( l ), which involves only P j ( and P x . Two spe- 
cific reasons for the failure of (I ) and (3) in the COT current sheet 
were identified. At Izl < z o the magnetic field direction was found to 
be irrelevant to the orientation of the calculated pressure tensor. 
The pressure tensor became field aligned near Izl = 1.5 z Q . Beyond 
here a small region was found in which the calculated P is field 
aligned but the pressure is dependent upon the azimuthal angle 
about the magnetic field direction. Equations (1) and (3) could be 
used only beyond Izl = 2z 0 where P is both field aligned and sym- 
metric about B. 

A thick quiet time model field was used for one COT study. This 
model was selected so that it could to be generated by an isotropic 
distribution function. As predicted by Wolf and Pontius [1993] and 
Usadi et al. (1996], currents were well approximated by (1) in the 
isotropic current sheet. The distribution functions of each of the 
separate ion groups, which were combined to produce the COT 
model, were anisotropic. As a result, (1) did not produce a good 
approximation to the current carried by individual ion groups, 
while (2) yielded the correct j y (*, z) for both individual groups 
and for the full consistent plasma. A thin current sheet model also 
was introduced. It was possible to generate a consistent thin current 
sheet with an anisotropic distribution function using only Speiser 
particles everywhere except in a region with adiabaticity parameter 
K near 0.8. 

The above nonguiding center results are important because 
highly structured distribution functions have been observed in the 
magnetotail [Nakamura et a!., 1991; Angelopoulos et al., 1992; 
Frank et al., 1996]. Our results suggest that one must evaluate the 
full pressure tensor to estimate forces in the inner current sheet. 


4.2. Substorm Effects 

Results from the COT model can be compared to previous anal- 
yses of the magnetotail. The Burkhart et al. [1992] 1-D simulation 
is particularly interesting. As in any 1-D study, the dP X2 /dz term 
generated the x component of particle forces, and the dP zz /dz 
term generated the z component. The ion source was located far 
from z = 0, so all particles started on Speiser orbits. As in our case, 
the off-diagonal pressure tensor element P xz was an order of mag- 
nitude smaller than the diagonal element P zz . Forces were well 
balanced in thin current sheets with small k. 

Burkhart et al. [1992] concluded that solutions could not be 
found in 1-D simulations when the average k of injected Speiser 
orbit ions approached one. Problems were encountered balancing 
the z components of particle and field forces. The COT method 
found solutions with k near one only in thick 2-D current sheets 
which contained many ions on trapped and cucumber orbits. 
Burkhart et al. [1992] noted that the current sheet tended to 
become thicker as k increased from small values even though all 
ions were injected on Speiser orbits. They concluded that there will 
be a catastrophic loss of equilibrium if K slowly increases in a thin 
current sheet. This is the same effect that is illustrated in Figure 6, 
where we were unable to find particles with k near 0.8 that could 
carry current in a thin sheet. The COT analysis found that other 
problems can arise in a thinning current sheet which is dominated 
by resonant particles. 

The following sequence of events, suggested by the above find- 
ings, may explain why many substorm onsets take place so close to 
the Earth, near midnight, and in a thin current sheet. A certain lobe 
field is needed so that the outward force exerted by the lobe on the 
magnetopause will balance the inward force exerted by the solar 
wind. A certain electric current is needed in the plasma sheet and 
magnetopause to generate the required lobe magnetic field. If the 
current sheet strength decreased in a small region of the tail, it is 
possible that more current could flow in adjacent regions to gener- 
ate the lobe field. However, if there are not enough particles or if 
the particles drift too slowly in a large region of the tail, something 
more dramatic must happen. A substorm onset which involves the 
collapse of a segment of the tail is one such dramatic event. Tail 
collapse will bring more ions into the region with insufficient 
cross-tail current, will change the thickness of portions of the cur- 
rent sheet, and will change k and the associated average cross-tail 
drift speed in the region of interest. The magnetopause shape and 
therefore the normal force exerted by the solar wind also will 
change. Another possibility is that guiding center electrons may 
become the dominant cross-tail current carriers in an extremely 
thin current sheet [Pritchett and Coroniti, 1995; Ma et al., 1995; 
Birn et al., 1996]. 

To be more specific, consider the T89 magnetotail which is a 
widely used model containing a thick current sheet. Assuming that 
5-keV protons carry the current in the T89 model, K is smaller than 
0.8 beyond about 12 R E at midnight, and k increases away from 
midnight [Pulkkinen et al., 1992; Kaufmann et al., 1993]. The 
cross-tail drift speed V. was found to be strongly dependent on the 
adiabaticity parameter k only in thin current sheets which are dom- 
inated by particles on Speiser orbits. Our COT analysis found that 
it is possible to generate thick steady current sheets with k = 0.8. 
The K parameter decreases as the tail becomes more stretched, B z 
decreases, and the current sheet becomes thin. Burkhart et al. 
[1992] and our COT analysis were unable to create a thin steady 
current sheet with K = 0.8. 



KAUFMANN ET AL.: NONGUIDING CENTER MOTION AND SUBSTORM EFFECTS 


22,167 


A substorm therefore would start at whatever location and time 
the current sheet thickness decreases to L 2 = z Q and k of the princi- 
pal ions also approaches 0.8. This point will be at midnight in T89 
because this is the local time at which the current sheet is thinnest 
and k is smallest. The onset location will be earthward of 1 2 /? if 
the T89 model correctly describes the magnetotail before thinning 
begins because K must initially be laiger than 0.8 if it is to decrease 
to 0.8. Figure 6 also shows that if the current sheet becomes thin 
cross-tail drift velocities become small over quite an extended 
region of the tail, centered at the point k = 0.8. Mitchell et al 
[1990] observed one current sheet that became thin and was domi- 
nated by Speiser particles just before a local disruption. Lui et al. 
[1992] found a series of disruptions in data from the AMPTE/CCE 
satellite which had an apogee of 8.8 R £ . The disruptions took place 
in current sheets with estimated thicknesses of 0.2 to 0.5 R and in 
regions with estimated adiabaticity parameters of k = 0.5 t<fo.9. 

Some mechanism must reduce the plasma content of the tail as 
it slowly becomes thin [Hones et al., 1971], It is generally thought 
that changes in the electric field or equivalently in the convection 
pattern are most important during substorm growth phase. A sud- 
den increase in the normal force exerted by the solar wind on the 
magnetopause could produce a small decrease in plasma sheet 
thickness. The pressure can respond both to changes in the solar 
wind momentum density and also to a wide variety of changes in 
the IMF direction. These small pressure changes could be impor- 
tant because a pressure pulse propagates rapidly to the inner mag- 
netotail near midnight where onsets appear to begin [Lopez et al., 
1990; Lui, et al., 1992], Lyons [ 1996] emphasized how frequently a 
change in the IMF and the associated convection pattern is associ- 
ated with substorm onset. The suddenness and localization of sub- 
storm onsets and the observation of multiple expansions imply that 
even triggered substorms involve some internal instability which 
produces one or a sequence of localized disruptions in the cross-tail 
current and its diversion to the ionosphere. 

Kaufmann et al. [this issue] stressed the possible significance of 
a net tailward force and the resulting stretching of the plasma sheet 
everywhere outside a substorm current diversion loop immediately 
after onset. This process will decrease both k and the current sheet 
thickness on the dawnside and duskside of the original diversion 
loop. Since k and the current sheet thickness tend to be larger away 
from midnight, these parameters must both decrease more at the 
dawn and dusk locations than at midnight in order to reach the L ~ 
z o and k = 0.8 values. The stretching after an onset therefore may 
cause a disruption region that starts near midnight to expand both 
toward the dawn and dusk. The result could be either a smoothly 
spreading or a multiple expansion substorm. 

The other distinct nonguiding center effect that was discussed 
also requires a thin current sheet dominated by Speiser orbits The 
average nonguiding center particle is already moving along the y 
axis at nearly its maximum speed in a thin current sheet dominated 
by Speiser particles with a k: near one of the resonances, such as k 
- 0.5. Unlike guiding center particles, the cross-tail velocity of a 
resonant particle cannot increase significantly when the current 
sheet gets still thinner and field lines become more sharply bent. 

The result again is an inability for ions to carry the current needed 
to produce the lobe field that must be present to balance the inward 
force exerted by the solar wind. 

The k = 0.8 chaotic particle and the ic = 0.5 resonant particle 
mechanisms described above both involve a thin current sheet. In 
both cases the existing ions may be unable to carry the necessary 
cross-tail current. The principal distinction is that disruption with k 


= is based on a drop in the average cross-tail velocity. Although 
the ion content of a section of the tail decreases as the tail becomes 
thinner [Hones et al., 1971], the K ~ 0.X mechanism is not primarily 
dependent on a decrease in the total number of ions. In contrast, the 
K = 0.5 mechanism is based upon a loss of plasma from the onset 
region. 

One other mechanism involving electrons which does not 
require a sudden change in plasma content also was discussed and 
could contribute to the disruption of cross-tail current. The precipi- 
tation of energetic electrons and their replacement by low-energy 
ionospheric electrons to maintain charge neutrality provides a 
mechanism that can reduce the electron component of the cross-tail 
current. The original energetic electrons drift faster and carry more 
current than the low-energy replacement electrons. 

The COT calculations assume steady state conditions, and those 
presented here were two-dimensional. A qualitative discussion of 
ways in which these results could fit into a 3-D time dependent 
substorm model was included. This was primarily a time depen- 
dent extension of the current diversion process discussed by Kauf- 
mann et al. [1994], The specific example considered was an 
evolution from one steady state condition to another. It was seen 
that ion and electron currents each are continuous in a steady state 
model. Electrons carry small steady field-aligned currents at the 
edges of a region of high or low density that is confined in the 
cross-tail direction. These currents are closed by cross-tail electron 
drift near the equator. Electron and ion currents are not individually 
continuous during the transient period in which the plasma density 
is increasing or decreasing. This situation requires equal net elec- 
tron and ion currents to flow into or out of the region of changing 
density [Atkinson, 1984], The result is a large field-aligned electron 
current which provides continuity during the transient phase when 
combined with the discontinuous ion cross-tail current. 

Many other mechanisms could produce substorm onset in a thin 
current sheet but were not discussed because they have been 
reviewed recently [Baker et al., 1996; Birn et al., 1996; Lui, 1996; 
Lyons, 1996; Rostoker, 1996; Sergeev et al., 1996]. The ’nonguiding 
center mechanisms discussed here are speculative. However, they 
require a specific current sheet thickness of about 2z 0 , an anisotro- 
pic distribution function dominated by Speiser particles, and either 
K near 0.8 or k near a resonance (e.g., k = 0.5) with large V /v 
and so should be testable. y 0 
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